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Abstract 

The Gaussian Wave-Packet phase-space representation is used to show that 
the expansion in powers of fi of the quantum Liouvihe propagator leads, in 
the zeroth order term, to results close to those obtained in the statistical 
quasiclassical method of Lee and Scully in the Weyl-Wigner picture. It is 
also verified that propagating the Wigner distribution along the classical tra- 
jectories the amount of error is less than that coming from propagating the 

Gaussian distribution along classical trajectories. 
PACS numbers: 34.10.+X, 03.65. Sq 
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I. INTRODUCTION 

Nowadays, the advantages and difficulties for using the phase-space formulation of quan- 
tum mechanics are well known ll|-^. This formulation remains still very useful for studying 
the classical limit of quantum mechanics as well as for describing semiclassical approx- 
imations in collisional processes; for these purposes many authors use the Weyl-Wigner 
(WW)picture |jl],^-0. For the collisional problem it is often combined with the quasiclassi- 
cal method of Ref. 0. Following this approach Lee and Scully |l|] have improved the accuracy 
of this method successfully with their Statistical Quasiclassical (SQC) method which was 
first suggested by Heller . As we have shown in a previous paper , the approach of Lee 
and Scully corresponds to the zeroth order term of the expansion of the quantum Liouvillian 
in powers of h in the WW picture. So, their calculated transition probabilities could have 
higher order corrections. 

The aim of this article is to show that the use of the Gaussian Wave Packet (GWP) 
phase-space representation |^,|TU|,|TT| gives for the zeroth order term, which corresponds to 
what we call the Causal Approximation (CA), results similar to those obtained in the WW 
representation '^. 

In the derivation of the expansion of the quantum Liouvillian in the GWP picture all 
orders of h are included; the ffist order term we shall call the Quasicausal Approximation 
(QCA). As it is shown in Refs. [|l3)0] one of the advantages of the GWP representation is 
that the quantum fluctuations of the mapped physical quantities become more evident; also 
the distribution function obtained when the density operator for a pure state is mapped into 
this representation is always non-negative. 

Here we are also interested in verifying numerically, for the collisional problem, the state- 
ment put forward by Lee 0: "the amount of error arising from propagating the Wigner 
distribution function (WDF) along the classical trajectories is usually considerably less than 
that coming from propagating other distributions along classical trajectories" . In our case 
the other distribution is the GWP distribution ||To|-[l2|. It is in agreement with Lee's state- 



ment the result we derive in Sec. II: first order corrections in h in the SQC method do not 
improve the numerical results of Lee and Scully |I| . 

In our comparison of the Weyl-Wigner formalism with the GWP approach we are going 
to work with a colhnear non-reactive collision of an atom with a diatomic molecule, the 
interaction between them being an exponential repulsion, first used by Secrest and John- 
son 



|13| , of the form V = Voexp[— a(x — y)] where the constant Vq is related to the classical 
turning point of the trajectory of the particle but its value does not have any effect on the 
results. In the WW and GWP formalisms we have taken this potential, the definition of 
the X and y coordinates, as well as the value of a and Vq from Ref. [|13|, the parameter a 



is adjusted by published experimental data [|T^]. The mapped Hamiltonian in the GWP 
formalism has the constant Vq renormalized. 

In Sec. II we shall summarize our previous paper (Ref. |Q). In Sec. Ill we introduce 
the GWP representation and the CA and QCA. In Sec. IV we derive expressions for the 
transition probabilities and numerical results are presented in Sec. V. 

II. CAUSAL AND QUASICAUSAL APPROXIMATIONS IN THE WW 

FORMALISM 

In the first part of this section we are going to give a review of Ref. [^. The quantum 



Liouville equation in the WW picture is |T^ 

dW{q,p,t) 

— = -iCQW{q,p,t), (1) 

where {q,p) is a point of phase-space and W{q,p,t) is the Wigner distribution function 
(WDF). We are using just one dimension. The quantum Liouvillian is 
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^Q = H{q,p) 



H{q,p) being the Hamiltonian of the system and the operator 



(2) 



^^ d_d_ d_d_ ,„x 

dq dp dp dq 



is the Poisson bracket, arrows indicate on which side the derivatives operate. 
In Ref. we show that from the formal solution of Eq. (0) given by 

W^(g,p,t)=e-*^«(*-*»)W^o(g,p), (4) 

where VFo(g,p) is the WDF at the initial time to- Taking the classical limit of Cq in Eq.(^ 
we get 

Cd = iH{q,p) A, (5) 

and Eq. (|) becomes 

W^''\q,p,t) = e-'^^'^'-''^Wo{q,p) = Wo{q{to-t),p{to-t)). (6) 

Thus, each point (g, p) of the phase space of the initial WDF evolves classically according 
to Hamilton's equations, following a classical trajectory reversed in time. This we call the 
Causal Approximation (CA). 

Still, according to Ref. |^, we can make an expansion of Cq in a power series of h'^ which 
is substituted into Eq. (P giving 
dW{q,p,t) 



+ t/:oW{q,p,t) = -tY.h'''£2nW{q,p,t), (7) 



dt 



where Cq = Cd and 



C2n = H{q,p) 



(-1)" /^\2„,+l 

22"(2n + l)! \^) 



(8) 



The integral equation corresponding to Eq.(^ is given by 

OO „f 

W{q,p,t) = e-^^»(*-*°)iyo(g,p) - ^ E ^'" / rft'e-*^°(*-*')/:2nVr(g,p,0- (9) 



n=l -'^0 

Solving this equation iteratively, we get in first order the Quasicausal Approximation (QCA) 

ft 

'to 

In Eq.(|TUp the operator exp[— ■j£o('^ — io)] is responsible for the classical character of the 
evolution between different times. 



WQCA{q.P.t) = e-'^°^'-'°'>Wo{q,p) - ih^ f' dt'e~'^''^'-''^C2e-'^°'^''-'°^Wo{q,p). (10) 

Jtn 



This formalism is applied in Ref. to a collisional process where a molecule suffers a 
collision from a pointlike projectile. As a result the molecule is transferred from the initial 
discrete energy level \i) to the final level |/), the total probability for this transition, in the 
limit t ^ oo and to -^ — oo, is given by 



/oo roc 

-oo J ~oo 



(11) 



Eq.(|TT]) is an exact result. Introducing now the QCA we get in the limit t ^ oo and 



n 



QCA 



27lh 



dpo / dqoWi{qo,po)Wf{q{t - to),p{t - to)) 



dpo / dqoW,iqo,Po) / rftV^^^*"* ^C2Wfiqit' - to), pit' - to)) 

-oo J —oo Jto 



(12) 



where the first term corresponds to the CA, while the second one is the QCA. The CA 
corresponds to the statistical quasiclassical (SQC) method of Lee and Scully given in Ref. [|I|. 
Here as in Ref. [|I| the H2 — He collision is considered, H2 and He being treated as an 
harmonic oscillator and a free particle respectively. The Hamiltonian in the Weyl-Wigner 
phase space is given by 



H{Q,q,P,p) 



2M 2m 



(13) 



where Q and q are the translational and vibrational coordinates respectively, P and p being 
their respective momenta. All the parameters appearing in Eg. (p!3D , M,m,VQ,a and k = 



rauo 



(elastic constant of the oscillator) are taken from Ref. ||13| 



Following Lee and Scully (Refs. [||,|1) Q,q,P,p obey Hamilton's equations, so they de- 
scribe classical trajectories, the initial state Wi{qo,po) in Eg. ([T2|) , is given by the WDF for 
the harmonic oscillator, the pair (go, Po) refers to the initial position and momentum of the 
harmonic oscillator and it belongs to a two-dimensional rectangular grid whose size and 
density depend on the desired accuracy. 

Integrating numerically Hamilton's equations for each (go, Po) of the grid for the harmonic 
oscillator and the appropriate initial Q and P of the particle, we get the set of final pairs 



{q,p) for the oscillator and final {Q,P) for the particle, initial and final Q must be taken 
sufficiently large so that the particle can be considered free, which can be verified by using 
the fact that the total energy must be conserved along the trajectories. 

In Eq.(|l^) the final state of the system is now given by the WDF, Wf{q,p) calculated 
for all final phase space points of the grid. Once initial and final WDF are calculated for 
each point, the CA can be obtained using the first term in Eg . ([12|) . 

In order to obtain the QCA, the second term in Eq.(|l^) must be calculated. This term 
can be approximated by |]12| 

'to 
with the constant B = 7r(/ia)^Vo/12 and 



Cqca = BF,f f rft'e-"[«(*'-*«)-'?(*'-*°^ (14) 

Jtn 



dp, / rfgoW^i(go,Po)e'^*-*«)^«^T^iy/(go,Po). (15) 

-OO J —OO ^Po 

Making y = 2EoH/ft-^ {Eoh{(1o,Po) being the energy of the classical harmonic oscillator) 
we have 

g-^ Wfiqo, Po) = AopI + Bopo, (16) 



where 



xnuomj oy-^ xnuomj oy^ 



being w/(y) = lV/(go,Po)- 

Now we are going to show that this correction, given by Eq.(|l4|) oscillates periodically 
in the time t. 

Let Im be the time at which the atom is considered to be a free particle after colliding 
with the molecule. At this time the coordinate and momentum of the classical harmonic 
oscillator are given by 

qoM = A cos{ujtM + 0o) and poM = -mujA sin(co'tj\/ + 0o)- 

At a later time t = Im + At, with At > 0, the momentum po(^) of the harmonic oscillator 
will be 



Po{i) = Pom cos(c<j At) — mcuqoM sm{uj At). (17) 

Thus, making a = poA/ and b = —mujqoM, one has 

pl(t) = a^ cos^ LuAt + b^ sin^ iuAt + Sab"^ cos ujAt sin^ uAt + Sa^b sin cuAt cos^ a; At. (18) 
From Eq.(|l3) and Eq.(ll|) we obtain for Eq.(|l6D 

^oPo(^) + -BoPo(^) = 03 cosScuAt + fessinScuAt + Oi coscuAt + 61 sincuAt, (19) 

where 

as = Aq (a^ + a*^) , b^ = iAo{a^ - a*^) , 

Oi = (SAoaa* + BQ){a + a*), 61 = i{3AQaa* + -Bo)(tt — a*), 

and a = (l/2)(a — ib). Eq.([T9|) is a Fourier's series, which substituted into Eq.(|15D shows 
that Cqca in Eq. ([T^ ) is a periodical function, since in this equation the integral in the time 
converges in the limit t -^ 00. The average of this periodical function over one period of the 
oscillator will be zero. In this derivation we have used the approximate expression for the 
correction of the transition rate given by Eq.(|I^). The exact demonstration, although more 
envolved, follows along similar lines. 

III. CAUSAL AND QUASICAUSAL APPROXIMATIONS IN THE GWP 

FORMALISM 

In the GWP representation [|T^-|T2[, operators can be mapped both into a covariant {CV) 
and a contravariant (CTV) form and there are expressions which relate the CV with the 
CTV form as well as both of them with the corresponding WW representation. 

The commutator of two operators of the Hilbert space A and B in the CV form is 
written ||,[TD|jn| 

{pq\[A, B]\pq) = ^^^(g,p) (v - f) B'^''{q,p), (20) 



where (g,p) is a point in a phase-space, \pq) represents the minimum uncertainty gaussian 
wave-packet or coherent state, A^^{q,p) and B'"^{q,p) are the CV forms of operators A 
and B and r= exp[{h/2) D -D*], where D = {l/ao)d/dq — iaod/dp. Arrows indicate on 
which side operators act and qq is a constant with dimensions M^T'^, 

The quantum Liouville equation in the CV representation for the density operator 
P^^{q,p,t) = {pq\'^it)){^{t)\pq) is given by 

d 



dt 



P^'{q,p,t) = -^C^^P^■-{q,p,t) 



C being the quantum Liouvilhan in the CV form 



£^^ = ^^^^(g,p)(r-r 



and 'MP^{q,p) = {pq\H\pq) the CV form of the hamihonian H [|rT[. 

Now expanding C^^ in a ^ power series, C'~^^ = S^o ^"-^n; "^^ identify 

n+r 



C-. 



2"(n + l)! 



Im 



J^n+l^CV 



D* 



where Cq = Cd is the classical Liouvilhan. 
One defines a Green's function []TT| by 



M;(g,p,t|go,Po,^o) = e '^* *"^^ w(g,]?, to|go,Po, ^o), 



(21) 



(22) 



(23) 



(24) 



with the condition 



limw(g,p,t|go,Po,^o) = ^{q - qo)^{p - Po)- 

t^to 



This Green's function satisfies the Liouville equation 
d 



dt 



wiq,P,t\qo,Po,to) = -iC w{q,p,t\qo,po,to), t > to 



(25) 



(26) 



where the pairs (g^Po) , and {q,p) are the momenta and the coordinates at times to and t 
respectively. Here, as discussed in Ref . , classical causality is broken in the h power series 
expansion of C'"^ when terms with n > 1 are retained. 
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Because of Eq. (|5), Eq. ( p6D can be written as follows 



Q 2 °° fh\n+l 



(27) 



In the GWP phase space the CV representation of the density operator, P'-^^(q,p,t) is 
given by ini 



/OO rOQ 

dpo / dqow{q,p,t\qo,Po,to)P'^^{qo,Po,to), 
-OO J —OO 



(28) 



so, if a formal solution of Eq. (p 



w{t\to) = e~'^(*-*")^"w(to|to) + 



9 rt OO I ii I 



(29) 



is substituted into Eq. (p8|), keeping in mind that u;(to|io) = ^(^ ~ Q'o)^(p ~ Po), the first 
term gives for P'-''^{q,p,t) the causal approximation (CA), while by including the lowest 
correction n = 1 term, the quasicausal approximation (QCA) is obtained. 

IV. TRANSITION PROBABILITIES FOR A COLLISIONAL PROCESS 



When we have a system in a given initial state \i), at time to and final state |/) at time 
t, with t > to, the transition probability Vi^j may be written [IT 



\ /*oo /*oo 

ZTTit J —OO J —OO 



ZTX I i J —OO J —OO 

]^ /"OO /"CXD 

27r^ 



dpi dqPr''iq,P,to)Pf''iq,P,t) (30) 

z,/l /4 J ~oo J —OO 

being Pf^^ the CTV distribution function of the initial state and P9^ the CV distribution 
function of the final state. 

For the CTV distribution function we have an equation |]ll| analogous to Eq. (|2lD whose 
formal solution is given by 



CTV 



Pr'{q,P,t)=exp -.(t-to)£^^^ Pr^'iq.P.to 



yCTVi 



10 



with 



and 



h 



^CTV /^CTV\ *■ 

r - r 



r =expi-h D D /2j . 

We use the model and method described by Lee and Scully |l],|| for the one-dimensional 
atom-molecule coUinear collision, but in the GWP phase-space formulation. Here also in 
this non-reactive process the atom is treated like a free classical structureless particle while 
the molecule is represented by a harmonic oscillator. The interaction between the atom and 
the molecule [^ is the exponential repulsion described in Sec. I. The Hamiltonian for this 



system in the CV representation is 

w"' = ^ + = + -M + v.//^-'"-" + 7 1 5: + 3 I . (31) 






where Q and q are the translational and vibrational coordinates respectively, P and p their 
respective momenta, K// = Voexp^a'^h/Aal) and the parameters M, m, Vq, a and k (elas- 
tic constant) are taken from Ref. |T3|, fitting the He-H2 system. The CTV form of the 
Hamiltonian is obtained from Eq. (^) by replacing Oq by —a^ |Tl|| . 
In the h power-series expansion of C^'^ 

Cpy = £0 + ^1^1 + ^^1^2 + ■ ■ ■ (32) 

the contributions of the coordinates Q and P were neglected, except in the zeroth order 
term given by 

°~'l^ dQ dP dP dQ^ dq d^~ dp d'qj ^ ' 

and which corresponds to classical motion. 

In this approximation D'^'^^Ti^^ = {a / ao)"'~^^Vef fexp[—a{Q — q)] for n > 1, if one takes 
into account that a^ = muj. Eq. ( p9|) can then be written 



Ti Jto 



E ^^^V^^(i'l^o) 

(34) 
11 



Substituting Eq. (p4D into Eq. (^) we have in first iteration 



^''''(g,P,t) = e-^(*-*°)^°P^^(g,p,to) 



+ 



2K 



-^ / dt' 



^-i(t-t')Co^-a{Q^q) 



ha T~)*\n 






n=2 



In order to obtain P'-^^{q,p,t) = exp[i(t — tQ)C'^^]Pj^^{q,p,to), which is the function 
appearing in Eq. (^0]), we proceed in the same fashion, but now as the Green's function is 
given by w{t\to) = exp[i(t — to)C'-''^]w(to\to), we get instead of Eq. ( ^51) 

P^^(g,p,t) = e^(*-*«)^«P^^(g,p,to) 



^'" ^ dt' 



h Jto 

since P^^{q,p,to) = P^^{q,p,to 



^i{t-t')Co^~a(Q~q) 



oo ( ha r)*\n 
n=2 ^ 



n\ 



Substituting Eq. (|36|) into Eq. (^) we have, with quantum corrections in all orders of h 
dpdq 



V,^, = l^lf^'-'(g,p,t,)e'"-"-"^'Pf"(g,p,t„ 



aQ J ZttTi Jto 



g«(i-i')-Cog-a(Q-g) 



• e*(*-*')^° Im 







'd^e^^' - l) D*e^(*'-*o)^"P.^^(g,p,to) 



(37) 



where we have used the property, valid for any operator A 

rl °° An 

Jo ^An + l)\ 



(38) 



n=0 



The first term in Eq. (^) is the CA which is the zeroth order term of our h power series- 
expansion, and corresponds to the expression of the Statistical Quasiclassical (SQC) method 
of Lee and Scully ||1|. The difference between these two expressions lies in the distribution 
functions, while Lee and Scully work with the product of two Wigner Distribution Functions 
(WDF), in Eq. (^) we have the product of two Gaussian Distribution Functions (GDF), 
one of them in the CTV form and the other in the CV form. 

Like in the SQC method, |I|J^ here Pf^^{q,p,to) represents the initial vibrational state 
of the system, which we propagate along the classical trajectories. We also use the Lee and 
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Scully method for constructing the two-dimensional rectangular grid in the {q,p) plane. 
Pf^^{qn,Pn, to = — C)o) is the weight carried by the point {qn,Pn) of the nth cell of the grid. 

After integrating Hamilton's equations for each point of the grid, the final GDF 
{P9'^^{q,p,t = oo)) is calculated and then the transition probabilities are computed in 
the CA which corresponds to the first term in Eq. (|37D. 

The GDF in the CV form is given by (see Appendix) 

P^^''iy) = ^re-y, (39) 

where y = E/hw, E{q,p) being the energy of the classical harmonic oscillator, n is the quan- 
tum number which corresponds to the nth eigenstate of the harmonic oscillator. Introducing 
the new variable r^ = y, the GDF in the CTV form is (see Appendix) 

pf^^(r)=e-^^^^Pf^(r). (40) 



V. NUMERICAL RESULTS AND CONCLUSIONS 

Using the mapped form of the quantum Liouville equation into the GWP phase-space 
we have derived expressions for the transition probabilities for semiclassical calculations of 
inelastic atom-molecule collisions which include also the quantum terms. These expressions 
are similar to those obtained in Ref. |^ for the WW formalism. 

We show in Tables I and II the transition probabilities Vi^f from the initial (i) to the 
final (/) state of the He-H2 system which were computed first by integrating numerically 
the Schrodinger equation, (these are the exact quantum mechanical results (QM) taken from 
Ref. |13[)) second by using the Weyl-Wigner representation in the statistical quasiclassical 
(SQC) method taken from Ref. |I[] and third by using the gaussian wave-packet representa- 
tion in the causal approximation (GWP). The last ones are our results and they are given 
by the first term in Eq. (^7\). For more details on the system and models see the mentioned 
references. 

13 



The results in Tables I and II obtained using only the zeroth order term for both methods 
SQC and GWP should improve if quantum corrections were introduced. 

The accuracy of the SQC and the GWP methods was checked by studying how well 
the transition probabilities obey microscopic reversibility, the results are presented in Table 
III. These results confirm that for the CA, the amount of error coming from propagating 
the WDF along classical trajectories is less than that coming from propagating the GWP 
distribution also along the classical trajectories, as stated in Ref. 0. This agrees with the 
vanishing result found in Section II for the correction to the CA in the WW formalism. How- 
ever, preliminary calculations of the first quantum correction (QCA) in the GWP formalism 
gave nonzero results. 

Both methods (SQC and GWP) give nonvanishing results for classically forbidden pro- 
cesses. The numerical calculations, as far as the CA is concerned, for both formalisms 
present almost the same degree of difficulty for being performed. 

ACKNOWLEDGMENTS 

One of us (M. C. T.) would like to thank the Instituto de Fisica Teorica for the hospitality. 

APPENDIX: DERIVATION OF THE CTV FORM OF THE GDF 

Let |v?n), with n = 0, 1,2... be the eigenstates of the harmonic oscillator and \a) = \pq) 
a coherent state, the CV form for the density operator |p!0| , p!5| {pq\^n){^n\'P'i) is given by 



P^^{r) = ^e '' , where a = {al/2h)2q + i(l/2hal)2p = re*^. If i? represents the energy of 
the classical harmonic oscillator and r^ = -^ = y, we may write 

P^'^iy) = V"e-^. (Al) 

nl 

We see that this CV form of the GDF depends only on the energy of the harmonic oscillator. 
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In order to obtain the CTV form of the GDF we use the property of the coherent 



states |T^ which relates the normal {NO) and antinormal (AO) ordering of a function / of 
a and a* 



f(^o\a,a*) = e^f''''\a,a*). 



In polar coordinates r and 6, we have 

92 _ I ( d d 1 92 ^ 
dada* 4r \dr dr r dO"^ 

and Eq. ( |A^) in the GWP representation is written llO 



:>CTV, 



r,e) 



4r I Or Or r 



a-" 



-)Cy I 



For any positive integer n we have derived the formula 



4r dr dr I 



-irn!L„p)e- 



where L„ are the Laguerre polynomials. 

Since P^^{r) for the harmonic oscillator do not depend on 9, Eq. 

as 



p^^v^r)=e-^o^'e^Pyir). 



Using Eqs. (|A6|) , Eq. ( |Al|) and Eq. (|A5|) we obtain for n = 

00 
Pr''ir) = Y.Ln{r')e-''. 

n=0 

For n = 1 we get, similarly, 

00 

n=0 

because 
For n = 2 



yCTV 



P-^^(r) = -r^e- + 



E 

.n=3 



n(n — 1) 



L„(0-2En^n(r2) + E^^ 



(A2) 



(A3) 



(A4) 



(A5) 



can be written 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



n=2 



n=l 
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Table Captions 

Table I. Transition probability Vo^f for a collinear He- 7^2 collision calculated by the quan- 
tum mechanical (QM) method (Ref. ||13|)) statistical quasiclassical (SQC) method (Ref. [0) 
and gaussian wave-packet (GWP) method in the CA (first term in Eq.(^)). The total ini- 
tial energy E is measured in units of huj/2, where uj is the vibrational frequency of the H2 
molecule. In the calculation of P^^^ at least 100 Laguerre polynomials were used, although 
only about 10 polynomials are required in order to obtain convergence. Numbers inside 
brackets give the upper limit of the transition probability and * means that the transition 
is prohibited classically. 

Table II. Similar to Table 1; except that the oscillator goes from the first exited state to 
the final state /. 

Table III. Similar to Table 1; results of Table 1 and Table 11 are compared in order to test 
microscopic reversibility. 
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TABLES 



E 


Vo- 


-/ 


QM 


SQC 


GWP 


8 


0- 


>0 


(0.892) 


0.893 


0.830 




0- 


^ 1 


0.108 


0.107 


0.156 




0- 


>2 


0.001 


- 


0.014 


12 


0- 


>0 


(0.538) 


0.529 


0.501 




0- 


^ 1 


0.394 


0.412 


0.349 




0- 


>2 


0.068 


0.068 


0.122 




0- 


>3 


- 


- 


0.028 


16 


0- 


>0 


(0.204) 


0.187 


0.229 




0- 


> 1 


0.434 


0.422 


0.339 




0- 


>2 


0.291 


0.314 


0.250 




0- 


*3 


0.071 


0.077 


0.124 




0- 


>4 


- 


- 


0.045 




0- 


>5 


- 


- 


0.013 


20 


0- 


>0 


(0.060)* 


0.046* 


0.090* 




0- 


> 1 


0.128 


0.202 


0.221 




0- 


>2 


0.366 


0.351 


0.270 




0- 


>3 


0.267 


0.294 


0.220 




0- 


>4 


0.089 


0.106 


0.134 




0- 


>5 


- 


- 


0.66 



TABLE L 



E 


Vi- 


-/ 


QM 


SQC 


GWP 


8 


1 - 


>0 


0.108 


0.106 


0.135 




1 - 


> 1 


(0.850) 


0.863 


0.780 




1 - 


>2 


0.042 


0.031 


0.085 


12 


1 - 


>0 


0.394 


0.411 


0.396 




1 - 


> 1 


(0.244) 


0.176 


0.250 




1 - 


>2 


0.345 


0.385 


0.272 




1 - 


>3 


0.037* 


0.028* 


0.082* 


16 


1 - 


>0 


0.434 


0.420 


0.377 




1 - 


> 1 


(0.034) 


0.065 


0.137 




1 - 


>2 


0.220 


0.151 


0.176 




1 - 


>3 


0.261 


0.302 


0.181 




1 - 


>4 


0.051 


0.061 


0.098 




1 - 


>5 


- 


- 


0.031 


20 


1 - 


>0 


0.218 


0.199 


0.230 




1 - 


> 1 


(0.286) 


0.285 


0.233 




1 - 


>2 


0.009 


0.042 


0.143 




1 - 


>3 


0.170 


0.090 


0.149 




1 - 


>4 


0.240 


0.262 


0.149 




1 - 


>5 


0.077 


0.121 


0.105 



TABLE II. 
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E 


n- 


-/ 


QM 


SQC 


GWP 


8 


0- 


> 1 


0.108 


0.107 


0.156 




1 - 


>0 


0.108 


0.106 


0.135 


12 


0- 


> 1 


0.394 


0.412 


0.349 




1 - 


>0 


0.394 


0.411 


0.396 


16 


0- 


> 1 


0.424 


0.422 


0.339 




1 - 


>0 


0.434 


0.420 


0.377 


20 


0- 


> 1 


0.218 


0.202 


0.221 




1 - 


>0 


0.218 


0.199 


0.230 



TABLE III. 
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